Polarization multistability of cavity polaritons 
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New effects of polarization multistability and polarization hysteresis in a coherently driven po- 
lariton condensate in a semiconductor microcavity are predicted and theoretically analyzed. The 
multistability arises due to polarization-dependent polariton-polariton interactions and can be re- 
vealed in polarization resolved photoluminescence experiments. The pumping power required to 
observe this effect is of 4 orders of magnitude lower than the characteristic pumping power in 
conventional bistable optical systems. 

PACS numbers: 71.36.+C, 42.65.Pc, 42.55.Sa 



O 

o 



> 

00 
(N 

o 

-I— > 



T3 

c 

o 
o 



X 



Introduction. — Cavity polaritons are elementary exci- 
tations of semiconductor microcavities with extremely 



ranging from 10 to 10 of 



small effective mass m 
the free electron mass [l[. 

In the low density limit they behave as weakly interact- 
ing bosons and their Bose-Einstein condensation (BEC) 
has been recently claimed Q. Under resonant excita- 
tion two main nonlinear mechanisms have been identi- 
fied, the polariton parametric scattering 3, 4, [a, [6|, M, [a] 
and bistability of the polariton system [9|, [lCj, [ll|, [l2|, ll3| . 
These two nonlinear mechanisms often coexist, being of 
the same origin 11 , ■ 



Another important peculiarity of the cavity polaritons 
is related to their polarization or pseudospin degree of 
freedom. Polaritons have two possible spin projections 
on the structure growth axis, ±1, corresponding to the 
right and left circular polarizations of the counterpart 
photons. In case of non-zero in-plane wave- vector (k ^ 0) 
these two components are mixed by TE-TM splitting [3] , 
but we will be interested in k = case. The mixing of the 
k = states appears due to the polariton-polariton inter- 
action, which depends on the spin orientation. Namely, 
the interaction of polaritons in triplet configuration (par- 
allel spin projections on the structure growth axis) is 
different from that of polaritons in singlet configuration 
(antiparallel spin projections). The spin-dependent po- 
lariton coupling strongly affects the predicted superfluid 
properties of the polariton system [HI, [l6| and leads to 
remarkable nonlinear effects in polariton spin relaxation, 
such as self-induced Larmor precession and inversion of 
the linear polarization during the scattering act 



17] 



In this Letter we show that the interplay between the 
nonlinearity caused by the polariton-polariton interac- 
tions and the polarization dependence of these interac- 
tions results in a remarkable multistability of a driven 



polariton condensate, contrary to the usual optical bista- 
bility in the spin-less nonlinear case. We consider the 
ground state (k = 0) polariton mode having a finite life- 
time which is coherently pumped at normal incidence to 
the microcavity plane by a cw laser of variable intensity, 
frequency, and polarization. We calculate the resulting 
intensity and polarization of the polariton field and show 
that, for a given polarization of the pump and depend- 
ing on the history of the pumping process, the polariton 
polarization can, in general, take three different values. 
For instance, a linearly polarized laser can result in a 
strongly right-circularly, strongly left-circularly, or in a 
linearly polarized polariton state. The whole system is 
therefore not simply bistable but tri- or even multistable. 
It is worth noting that this conversion from linear to cir- 
cular polarization arises without any TE-TM splitting 
[HI and represents a new nonlinear effect. 

Qualitatively, the multistability can be understood as 
follows. Assume that the cavity is illuminated by a laser 
light at normal incidence at the frequency above the bot- 
tom of the lowest polariton branch (LPB). At low pump- 
ing, the pump is not in resonance with the polariton 
eigenstate, so that the population of the driven mode 
remains low. At higher pumping, polariton-polariton in- 
teractions lead to the blue-shift of the LPB, so that it 
approaches the pump laser frequency. At resonance, the 
population jumps up abruptly. If the pumping power 
is then decreased, the population of the polariton mode 
jumps down back, but at a lower threshold. As a result 
the typical S-shape dependence of the intra-cavity field 
on the pumping intensity appears as Fig. Q] shows, which 
means the formation of a hysteresis cycle. The additional 
polarization degree of freedom makes this picture much 
more complex and rich. At elliptically polarized pump- 
ing the hysteresises for right- and left-circularly polarized 
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FIG. 1: The condensate population \ip a \ versus pumping 
intensity \p<r\ 2 for the circularly polarized pump (a = ±1). 
Dashed line shows the unstable region. Ii and I2 mark the 
intensities of two jumps corresponding to the increase and de- 
crease of the pump, respectively. The numbers shown (scaled 
to the excitation area S ~ 100 /1m 2 ) are typical for cavity 
polariton experiments. 



components of the field do not coincide, which results in 
a sequence of polarization jumps and multistability. 

Coherent polarization evolution. — We consider a semi- 
conductor microcavity in the strong coupling regime 
pumped by a laser light perpendicular to the cavity plane 
(k = 0), with spatially uniform and slowly changing with 
time intensity p a and frequency uj near the bottom of 
the LPB, where a — ±1 is the circular polarization in- 
dex. We neglect the non-parabolicity of the polariton dis- 
persion, and the wave vector dependence of exciton and 
photon fractions. These approximations are well justi- 
fied since the effects we discuss are taking place in a nar- 
row energy and wave vector range close to the polariton 
ground state at k — with the eigenfrequency ojq. 

The k — polariton wave function i/j a satisfies the 
driven spinor Gross-Pitaevskii equation [ijj Il9j| . which 
in the quasi-stationary case can be written as (we set 

h= 1) 
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where r is the polariton lifetime, ai( 2 ) is the ma- 
trix element of polariton-polariton interaction in the 
triplet (singlet) configuration, respectively. The exchange 
interaction is suppressed for polaritons in the singlet con- 
figuration since it involves energetically split-off dark ex- 
citons as intermediate states. As a result, there is an at- 
traction between polaritons of opposite sign, a 2 < 0, and 
polariton-polariton interaction is strongly anisotropic 

\a 2 \ «ai m. 

At thermodynamic equilibrium the polariton BEC 
leads to the formation of the linearly polarized conden- 



sate, which minimizes the free energy of the system [15| . 
In the case of cw driven system the situation is qualita- 
tively different. The polariton system is intrinsically out 
of equilibrium and its state is now driven by the pump. 
The parameters of the driven mode, however, are not 
uniquely defined. E.g., with a linearly polarized pump 
the polariton condensate can change its polarization ei- 
ther to the left- or to the right-circular, increasing this 
way the blue-shift and entering into resonance with the 
pump. The left and right circularly polarized conden- 
sates can appear in this case with equal probabilities. 

Equation {T]) can be used to relate the occupation of 
the condensate N = 



iV'+il 2 + l^-il^ and the degree of 



its circular polarization p c = 
the pump intensity / = \p+i 
polarization degree p p = (\p+i 



(\1>+i\ 
2 + \p 



- \4>-i\ 2 )/N with 
-i| 2 and the pump 



\p-i\ 2 )/I. Namely, 



I /At = [n 2 + t- 2 + (01 - a 2 )(l - p 2 c )nN 

+ (l/4)(a 1 -a 2 ) 2 (l-p 2 c )N 2 ]N, (2) 

Pp I/At= [n 2 +r- 2 -{l/A){a 1 -a 2 ) 2 (l-p 2 c )N 2 ]p c N, (3) 

where Q = lo — loq — a\N . 

It follows from Eqs. ([2]) and ([3]) that only for the pump 
with a pure circular polarization, p p = ±1, the polar- 
izations of the polariton system and the pump coincide. 
In this case the present spinor model can be reduced to 
the scalar model 
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12| . For elliptical pumping the po- 
larizations of the polariton system and of the laser are 
different due to the different blue shifts for right and left 
circular polarized components. Due to the effect of the 
self-induced Larmor precession [13] the polarization el- 
lipse of the driven mode is rotated with respect to the 
polarization ellipse of the pump by some angle which de- 
pends on the pump intensity. Moreover, in this regime 
even the signs of circular polarization degrees of the po- 
laritons and the pump can be different, as we show in the 
next subsection. 

Multistability and hysteresis. — The polarization and 
occupation of the driven mode, p c and N, can be de- 
termined self-consistently from Eqs. ^ and ([3]), which 
can have several solutions depending on the value and 
polarization of the pump. 

The condensate population IV'iil 2 as a function the 
pump intensity |p±i| 2 for the case of fully circular po- 
larized excitation (when the polarization of the driven 
mode coincides with that of the pump) is shown in Fig. 
[U The energy of the laser u> is chosen above the energy 
of the bare polariton state, u) — ujq = 3 meV, so that the 
curve shows the classical S-shape. The matrix element 
cii = GxEba 2 ^/ S, where ob — 100 A is the two dimen- 
sional exciton Bohr radius, E\> = 8 meV is the exciton 
binding energy, x = 1/4 is the squared exciton fraction, 
and S = 100 /1m 2 is the laser spot area. The polariton 
life-time is r = 2 ps. These parameters are typical for 
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FIG. 2: (a) Circular polarization degree of the condensate 
versus external pumping intensity for slightly elliptical pump 
(p p = 0.2). Arrows show the direction in which the pump in- 
tensity is changed, (b) Circular polarization degree of the 
driven mode p c versus circular polarization degree of the 
pump pp. Arrows show the direction in which the pump polar- 
ization degree is changed. The pump intensity is just above 
It. Crosses mark the stability points, which correspond to 
p p = 0. 



a GaAlAs microcavity. We have denoted by 1\ the laser 
intensity corresponding to the bending point of the S- 
shape taking place with the increase of intensity, and 1% 
is the laser intensity corresponding to the bending point 
in the case of the intensity decrease. 

The evolution of the condensate polarization can be 
conveniently illustrated and understood considering the 
case «2 = 0, when Eqs. ((TJ for two values of a are sim- 
ply decoupled and the two circular components evolve 
independently 2l| . The change of the p c with the total 
pump intensity is shown in Fig. HJ^a) for the elliptically 
polarized pump with p p = 0.2, so that a + component 
slightly exceeds a~ one. As one increases the intensity 
of the pump, the solution moves along the lower branch 
of the S-curve for both a + and o~ components. How- 
ever, the a + component, since it dominates, reaches the 



threshold intensity I\ first. The corresponding intensity 
of the polariton field jumps abruptly to the upper branch. 
At the same time, the intensity of the a~ pump has not 
yet reached the l\ bending point. So, the jump of the 
total polariton density is accompanied by a jump of the 
circular polarization degree. If we increase further the 
intensity of the laser, the intensity of the a~ mode also 
reaches It. The polariton population increases again, 
but this now results in an abrupt decrease of the circu- 
lar polarization degree of the driven mode. If we now 
reduce the intensity, the reversed process takes place at 
the pumping intensity I2 so that hysteresises in both the 
occupation and polarization power dependencies appear. 

Fig. [2jb) shows another interesting configuration, 
where the laser intensity I is kept constant in the do- 
main I > I\ and It > I j1> 1%. The cyan line shows the 
change of the circular polarization degree of the driven 
mode p c induced by the laser initially polarized a + and 
whose polarization is progressively rotated towards the 
<7~ polarization. One can observe a weak decrease of p c , 
which however remains quite high even when the pump- 
ing is linearly polarized. This is due to the fact that the 
<7 + component remains on the upper branch of the S- 
curve whereas the a~ drops on the lower branch. Then 
there are two jumps of polarization corresponding to the 
jumps of cr + and o~ components of the polariton pop- 
ulation, and finally the polarization becomes fully o~ . 
The black line shows the evolution of p c with the in- 
verse change of the pump polarization. The stable points 
corresponding to full linear polarization of the laser are 
marked with crosses in this Figure. One can see that 
the internal polarization can be either nearly cr + , either 
nearly cr~, or fully linear. The latter case is in fact degen- 
erate. There can be two stable driven mode occupations 
N for the same value of the external laser intensity /. 

Fig. [3] shows the functional dependence between a + 
and a~ components of the polariton population and the 
intensity and polarisation of pump calculated accounting 
for the coupling between polaritons with opposite spins 
described by U2 — — O.lai. This value of 012 corresponds 
to recent estimations of Ref . [20] . The circular polariza- 
tion degree of the pump laser is represented by the color 
of the surface of solution and the intensity of the pump is 
on the vertical z-axis. The linear part of the polarization 
of the laser is kept aligned along x-direction. 

The green bright areas correspond to nearly linearly 
polarized pumping. If one increases the intensity of the 
pumping laser, keeping it linearly polarized, the system 
follows the black solid line and black arrows shown in 
the figure. One can see that from the critical point at the 
end of the solid black line the system can jump into three 
possible stable points (shown by crosses). One of them 
corresponds to linearly polarized state and two others 
to nearly right- and left-circularly polarized states. The 
choice of the final state by the system is random and is 
triggered by fluctuations. 
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FIG. 3: The pump intensity I and polarisation (color) versus 
the circular-polarized components of the driven mode IV'iil 2 - 
Color shows the polarization of the pump p p (bright green cor- 
responds to linear polarization). The crosses mark the four 
stable points for the driven mode corresponding to the same 
linearly polarized pump intensity. Arrows show the three pos- 
sible jumps in case if the pump intensity is slightly increased. 



Note also that the jump of the system to the circu- 
larly polarized final state induces a red shift of the cross- 
circular component because of the negative sign of ai. 
The red shift drives this component out of resonance with 
the pump, which leads to stronger polarisation of the final 
state. The positive feedback in polarisation of the con- 
densate would not take place for the case «2 = where 
only a linearly polarized driven mode would be formed. 
Experimentally, we expect this effect to have a key im- 
pact on the polarization measurements performed with 
resonantly excited microcavities. It will result in a ran- 
dom sign of the observed polariton polarization changing 
from one experiment to another. 

It should be noted that polarization multistability and 
chaos in nonlinear optical systems have been studied for 
more than 30 years (see, e.g., [22]). A number of in- 
triguing nonlinear effects have been predicted and ob- 
served in different systems, including four-wave mixing 
in anisotropic crystals 2^, 24 1. magnetic cavities 25 1 and 
VCSELs [26J . However, in the systems investigated pre- 
viously extremely high powers (10 MW/cm 2 [22j) are re- 
quired for observation of polarization multistability ef- 
fects linked to optical nonlinearities. The advantage of 
the microcavities is that in the strong coupling regime the 
nonlinear threshold corresponds to much lower pumping 
powers (e.g. 650 W/cm 2 



27] ) which really opens the 



pendent Gross-Pitaevskii equation. We have shown that 
for a given polarization of the pumping laser the polariton 
polarization can take three different values. For instance 
at linear pumping the resulting polariton state can be 
right circularly, left circularly or linearly polarized. In 
realistic cases, the driven mode polarization is random 
and is triggered by the fluctuations of the external laser. 

The multistability arises from the interplay between 
the anisotropy of the polariton-polariton interaction and 
the bistable behavior of cavity polaritons. 
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way to practical realization of new applications like data 
encryption and communication [28} . 

Conclusions. — We have analyzed the polarization of 
the spinor polariton system pumped by a cw laser of vari- 
able polarization at k = by solving the polarization de- 
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